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Abstract. Strong solutions of the non-stationary Navier-Stokes equations un- 
der non-linearized slip or leak boundary conditions are investigated. We show 
that the problems are formulated by a variational inequality of parabolic type, 
to which uniqueness is established. Using Galerkin's method and deriving a 
priori estimates, we prove global and local existence for 2D and 3D slip prob- 
lems respectively. For leak problems, under no-leak assumption at t = we 
prove local existence in 2D and 3D cases. Compatibility conditions for initial 
states play a significant role in the estimates. 



1. Introduction 

Let be a bounded smooth domain in M.'^ {d — 2,3), and fix T > 0. We 
suppose that the boundary T = dQ consists of two nonempty open subsets, that is, 
r = roUFi, Fonri = 0. We are concerned with the non-stationary incompressible 
Navier-Stokes equations in fl: 

u' + {u-W)u- jyAu + Wp = f in nx{0,T), (1.1) 

divu = in f7x(0,T), (1.2) 

with the initial condition 

u = uq in 17 X {0}. (1.3) 

Here, v, u, p, and / denote a viscosity constant, velocity field, pressure, and external 
force respectively; u' means the time derivative 

As for the boundary condition, we impose the adhesive b.c. on Tq: 

It = on To. (1.4) 

On the other hand, we consider one of the following nonlinear b.c. on Fi: 

U„ = 0, |crr|<5, ffr ■ Ur + g\ur\ — 0, On Fi, (1-5) 
which is called the slip boundary condition of friction type (SBCF), and 

Ur = 0, \crn\<g, CT„u„ +g|u„| = 0, on Fi, (1.6) 

which is called the leak boundary condition of friction type (LBCF). Here, n is an 
outer unit normal vector defined on F, and we write ii„ :— u ■ n and Ut ■— u — u„n. 
The stress tensor T = (T^ )i,j=i,...,d is given by = -pSij + i^(|fj + |^), % 
being Kronecker delta. We define the stress vector a = a{u,p) as a — Tn, and 
write cr„ := cr • n and ar ■= cr — Onn. One can easily see that cr„ — (7n{u,p) may 
depend on p, whereas CTt — ar{u) does not. 

The function g, given on Fi and assumed to be strictly positive, is called a mod- 
ulus of friction. Its physical meaning is the threshold of the tangential (resp. nor- 
mal) stress. In fact, if \ar\ < g (resp. |cr„| < g) then p.Sp (resp. p.6p ) implies 
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Ut — (resp. Un = 0), namely, no slip (resp. leak) occurs; otherwise non-trivial 
slip (resp. leak) can take place. We notice that if we make g = formally, (|1.5p 
and (jl.6l) reduce to the usual slip and leak b.c. respectively. In summary, SBCF 
and LBCF are non-linearized slip and leak b.c. obtained from introduction of some 
friction law on the stress. 

It should be also noted that the second and third conditions of (11.51) (resp. (HH])) 
are equivalently rewritten, with the notation of subdifferential, as 

CTr e -gdlurl (resp. o-„ e -gd\un\)- 

Though we will not pursue this matter further, one can refer to [51ll7| for the Navier- 
Stokes equations with general subdifferential b.c. See also [4], which considers the 
motion of a Bingham fluid under b.c. with nonlocal friction against slip. 

SBCF and LBCF are first introduced in [5J [S] for the stationary Stokes and 
Navier-Stokes equations, where existence and uniqueness of weak solutions are 
established. Generalized SBCF is considered in [121 HO]- The H^-H^ regular- 
ity for the Stokes equations is proved in j27j . In terms of numerical analysis, 
[a Ha m [21 EH El [25] deal with finite element methods for SBCF or LBCF. 
Applications of SBCF and LBCF to realistic problems, together with numerical 
simulations, are found in [T51 128) . 

For non-stationary cases, [7J[S] study the time-dependent Stokes equations with- 
out external forces under SBCF and LBCF, using a nonlinear semigroup theory. 
The solvability of nonlinear problems are discussed in |21j for SBCF, and in [1] for 
a variant of LBCF. They use the Stokes operator associated with the linear slip or 
leak b.c, and do not take into account a compatibility condition at t = 0. 

The purpose of this paper is to prove existence and uniqueness of a strong so- 
lution for (jl.ip - (ll.4p with (jl.5p or (II. 6p . We employ the class of solutions of La- 
dyzhenskaya type (see [18]), searching {u,p) such that 

u e L°°iO, T; H\n)'^), u' G L°°(0, T; L^{n)'^) f] L^{0, T; H'^ifl)''), 
peL°°(0,r;L2(r!)). 

There are several reasons we focus on this strong solution. First, from a view- 
point of numerical analysis, we would like to construct solutions in a class where 
uniqueness and regularity are assured also for 3D case. Second, we desire an L°°- 
estimate with respect to time for p, which may not be obtained for weak solutions 
of Leray-Hopf type (cf. [29, Proposition III. 1.1]). Third, in LBCF, it is not straight- 
forward to deduce a weak solution because of (II. 7p below. Similar difficulty already 
comes up in the linear leak b.c. (see [26) ') 

The rest of this paper is organized as follows. Basic symbols, notation, and 
function spaces are given in Section 2. 

In Section 3, we investigate the problem with SBCF. The weak formulation is 
given by a variational inequality, to which we prove uniqueness of solutions. To 
show existence, we consider a regularized problem, approximate it by Galerkin's 
method, and derive a priori estimates which allow us to pass on the limit to deduce 
the desired strong solution. Using the compatibility condition that uo must satisfy 
SBCF, we can adapt uq to the regularized problem, which makes an essential point 
in the estimate. 

Section 4 is devoted to a study of the problem with LBCF. There are two major 
differences from SBCF. First, as was pointed out in the stationary case ,B, Remark 
3.2], we cannot obtain the uniqueness of an additive constant for p if no leak occurs, 
namely, m„ = on Fi. Second, under LBCF, the quantity 

J ^^{u-V)v-v^dx^^ J u„\v\^ds (if divu^O) (1.7) 
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need not vanish because u„ can be non-zero. This fact affects our a priori estimates 
badly, and we can extract a solution only when the initial leak ||uon||L2(ri) is small 
enough. Incidentally, if we use the so-called Bernoulli pressure p+ instead of 
standard p, the mathematical difhculty arising from (ll.7|) are resolved; nevertheless 
the leak b.c. involving the Bernoulli pressure is known to cause an unphysical 
effect in numerical simulations (see [121 p.338]). Thereby we employ the usual 
formulation. 

Finally, in Section 5 we conclude this paper with some remarks on higher regu- 
larity. 

2. Preliminaries 

Throughout the present paper, the domain ft is supposed to be as smooth as 
required. For the precise regularity of fl which is sufficient to deduce our main 
theorems, see Remarks 13.51 and 14.41 We shall denote by C various generic positive 
constants depending only on ft, unless otherwise stated. When we need to specify 
dependence on a particular parameter, we write as C = C(/, uq), and so on. 

We use the Lebesgue space LP{il) {1 < p < oo), and the Sobolev space H^{n) = 
{(t)&L'^{n) I ||</'|ll^.(a)=E|a|<r <oo} for a nonnegative integer r, where 

if°(r2) means L^(Sl). H^{il.) is also defined for a non-integer s > (e.g. [lOl 
Definition 1.2]). We put i§(ri) ~ {q 'E L'^{Q) | J^qdx = 0}. For spaces of vector- 
valued functions, we write Lp(O)'^, and so on. 

The Lebesgue and Sobolev spaces on the boundary F, Fq, or Fi, are also used. 
i/°(Fi) means ^^(Fi), and we put //^(Fi) = {?/ e L^Ti) \ Jp^ ryds = 0}, where ds 
denotes the surface measure. For a positive function g on Fi , the weighted Lebesgue 
spaces Lg{Ti) and L!jY^(Fi) are defined by the norms 




respectively. The dual space of Lg(Fi) is L!jy^(Fi) (see Lemma 2.1]). 

The usual trace operator (jj i— >■ <^|r is defined from H^{D,) onto H^^^{T). The 
restrictions 4>\toj 4>\ti of (j)\r, are also considered, and we simply write (fi to indicate 
them when there is no fear of confusion. In particular, 7y„ and 1]^- means (7/-n)|r and 
(77— (?7-n)n)|r respectively, for 77 £ H^^'^{TY. Note that ||'7n|l_f/i/2(r) < C'||77||^i/2(p-)d 
and ll^^r ||i/i/2(r)ti ^ C'||77||^i/2(p)c! because n is smooth on F. 

The inner product of L^(il)'^ is simplified as (•,•); while other inner products 
and norms are written with clear subscripts, e.g., (•, •)^2(-pj-) or || • Wn'^in)''- For a 
Banach space X, we denote its dual space by X' and the dual product between X' 
and X by (•, ■) Moreover, we employ the standard notation of Bochner spaces 
such as L^{0,T;X), H^{0,T;X). 

For function spaces corresponding to a velocity and pressure, we introduce closed 
subspaces of H^(n)'^ or L^(f2) as follows: 

V = {ve H\n)'' 1 1; = on Fq}, V ^ {v e H^fl)'^ | u = on F}, 
Vn = {veV\ Vn = on Fi}, V^^ = {w e y I = on Fi}, 

To indicate a divergence-free space, we set H^{^Y — {v E H^{i}Y' \ div v ~ 0}. We 
use the notation ^ V n Hl{nY, = V r\ Hl{nY, Vn,a = K n H^i^Y, and 

Vr^„ = VrnH^^inY. 
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Let US define bilinear forms ao, h, and a trilinear form ai by 

ai(u,w,?«) = / {(-u • V)t;} • wda; (it, w, w G 

Jo. 

b{v,q) = - / divuqda; (u € H^{nf, q e 



The bilinear forms ao,& are continuous, and from Korn's inequality ([16| Lemma 
6.2]) there exists a constant a > such that 

ao{v,v)>a\\v\\l^f^,,^, (V« G F). (2.1) 

Concerning the trilinear term ai, we obtain the following two lemmas. 

Lemma 2.1. (i) When d — 2, for all u,v,w £ H^{il)'^ it holds that 

\ai{u,v,w)\ < C\\u\\]^i^^^^^\\u\\]^^^^^^\\v\\H^ (2.2) 

(ii) When d = 2 or d = 3, for all u,v,w G H^{n)'^ it holds that 

\ai{u,v,w)\ < C\\u\\]^^^^^y\\ufj^t^^y\\v\\Hii^^y\\w^^^ (2.3) 

Remark 2.1. In particular, we see from (12.31) that 

\ai{u,v,w)\ < C\\u\\Hi{n}AMm{ny\\M\m{<:iy- (2.4) 

Proof By the Sobolev embedding H^/^{n) C L^(fl) (resp. 77^/4(0) c L'^i^)) 
which is valid for d — 2 (resp. d — 2, 3), combined with an interpolation inequality 
between L^{il) and H^{n), we have 

ll"llL4(o)rf < C'||u||^i/2(n)d < C'll"lli2^(o)d||w|lHi%)<i = 2). 

(resp. ||u||L4(o)d < C||M||H3/4(n)d < C\\u\\]^2\i^-^A\u\\HHn)'^ (^ = 2,3).) 

Therefore, since |ai(M, u, i(;)| < C||u||^4(s-;-)d HwH^i^Q-jd ||t(;||^4(s-j-)d by Holder's inequal- 
ity, we conclude (resp. (gSl)). □ 

Lemma 2.2. (i) For a^/ m e y„,cr and w e H^{Q)'^, ai{u,v,v) — 0. 

(ii) For all u G yr.o- tf^c^ v G H^{il)'^ , ai{u, v,v) = ^ Jp^ w„|i;p c?s, and 

|ai(w,w,w)| < 7i||"«||L2(ri)||w||^i(s^)d, (2.5) 
where 71 is a constant depending only on fl. 
Proof. By integration by parts, we have 

ai{u,v,w) + ai{u,w,v) = — / divuv-wdx+ / UnV-wds, 



from which the conclusion of (i) and the first assertion of (ii) follow. Combining 
the Holder inequality \ai{u, v, v)\ < C\\un\\L^(ri) ll^l!L4(ri)<i' Sobolev embedding 
i?^/^(ri) C L'^iTi) {d = 2,3), and the continuity of the trace operator H^{fl) -^■ 
H^/^Ti), we derive ([231). □ 

Remark 2.2. Whether 71 is small or not, especially when compared to a in (|2.ip . is 
a very crucial point in our a priori estimates for LBCF (see Proposition [33J- This 
is why we distinguish 71 from other constants C and do not combine 71 with them. 
As (i) above shows, this problem does not happen when we consider SBCF. 



NAVIER-STOKES EQUATIONS UNDER FRICTIONAL BOUNDARY CONDITIONS 



5 



Furthermore, we introduce nonlinear functionals jr and jn by 

Jr{v)= f 9\v\ds (veL^r.f) and Mv) = f g\v\ds (7?eL2(ri)), 

where g > is a modulus of friction mentioned in Section [T] They are obviously 
nonnegative and positively homogeneous. In addition, they are Lipschitz continuous 
when g(t) e L^{Ti) for a.e. t g (0,r). 

The foUowings, which are readily obtainable consequences of standard trace and 
(solenoidal) extension theorems ([ini Theorems 1.1.5-6, Lemma 1.2.2], see also \W[ 
Section 5.3]), are frequently used in subsequent arguments. 

Lemma 2.3. (i) For v e Vn, it holds that ||'yr||Hi/2(ri)d < C'IK'llffi(n)''- 

(ii) For rj e H^^^(riY satisfying r]n = on Fi, there exists v € Vn,a such that 
Vr = ?/ on Fi and ||v||_f/i(n)'i < C||77||^i/2(ri)d. 

Lemma 2.4. (i) For v e Vr, it holds that ||t'n||iji/2(rj) < CH-wH jji(f2)d . 

(ii) For T\ e i/i/2(Fi) (resp. r\ G H^I'^{Vy)'^L\{Vx)), there exists v <^Vr (resp. v € 
K,cr) such that Vn = r\ on Fi and < C'||?7||^i/2(ri). 

The definition of a{u,p) given in Section [1] becomes ambiguous when {u,p) has 
only lower regularity, say u G H^{^lY, p G L^(il). Thus we propose a redefinition 
of it, based on the following Green formula: 

(— i/Au + Vp, t;) + J a{u,p) ■ V ds — ao{u,v) + b{v,p) (if divu = 0). 

Definition 2.1. Let u{t) e K, pit) € Q, u'{t) e L^{n)'^, f{t) e L^{ny. If ([LT]) 
holds in the distribution sense for a.e. t e (0, T), that is, 

{u\v) + aQ{u,v) + aiiu,u,v) + b{v,p) = {f,v) (Vu € V), (2.6) 

then we define a = a{u,p) e {H^^^iTiy)' by 

{cr,v)f^i/2(^j:'iy^^o{u,v)+b{v,p)~{F,v)y (Vw € F), (2.7) 

where F{t) £ V is given by {F, v)y — (/, v) — (u' , v) — ai{u, u, v). 

The above a is well-defined by virtue of the trace and extension theorem. It 
coincides with the previous definition when (u,p) is sufficiently smooth. In addition, 
by Lemmas US] and cr^ = cr-((T-n)n e {H^^^{Ti)'^y and cr„ = a-n e H^/'^{TiY 
are characterized by 

(CTr,Wr)Hi/2(ri)d = ao{u,v) + b{v,p) - {F,v)y^ (Vv E K), 
and 

(cr„,'y„)^i/2(r,) = ao{u,v) + b{v,p) - {F,v)y^ (Vw e K), 
respectively. By Lemma I2.3f ii) , ar actually does not depend on p. 

3. Navier-Stokes Problem with SBCF 

3.1. Weak formulations. Throughout this section, we assume f £L^{Qx {0,T))'^, 
uq G Vn^a, and g £ L^{Ti x (0,T)). Further regularity assumptions on these data 
will be given before Theorem l3.2l In addition, the barrier term jr is simply written 
as j. A primal weak formulation of (ll.ip ~ (ll.4p with (jl.5p is as follows: 

Problem PDE-SBCF. For a.e. t £ (0,r), find {u{t),p{t)) £ Vn x Q such that 
u'{t) £ L'^{VlY , u(0) = uq, ctt is well-defined in the sense of Definition 12. II \ar \ < g 
a.e. on Fi, and ar ■ Ur + g\ur \ = a.e. on Fi. 
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Remark 3.1. More precisely, "|(7r| < g" implies that ar G {H^^^{Ti)'^y actually 
belongs to L^^g{Ti)'^ with ||o-t-||^oo (p^^d < 1. In particular, ar € L^(ri)'^. 

Throughout this section, we refer to Problem PDE-SBCF just as Problem PDE. 
Similar abbreviation will be made for other problems. 

One can easily find that a classical solution of (|l.l|) - (jl.4[) with (|1.5|) solves Prob- 
lem PDE, and that a sufficiently smooth solution of Problem PDE is a classical 
solution. As the next theorem shows. Problem PDE is equivalent to the following 
variational inequality problem. 

Problem VI^-SBCF. For a.e. t E (0,T), find u{t) e Vn,a such that u'{t) G 
L^{n)'^, u{0) = uo, and 

(m', V — u) + ao{u, V ~u) + ai{u, U,V — u) +j{Vr) — j{Ur) >{f,V~ u) (Vw e Vn cr)- 

(3.1) 

Theorem 3.1. Problems PDE and VIo- are equivalent. 

Remark 3.2. The precise meaning of "equivalent" is that if solves Problem 

PDE, u solves Problem VIo-; if u solves Problem VIo-, there exists unique p such 
that solves Problem PDE. Hereafter we will frequently use the terminology 

"equivalent" in a similar sense. 

Proof. Let be a solution of Problem PDE. Then it follows that 

{u' ,v) + a(i{u,v) + ai{u,u,v) + b[v,p) - [ar,Vr)L^(TiY = {f^v) (Vw G Vn). (3.2) 
Using this equation together with |crT-| < g and CTt • Mt + 5|j^r| = 0, we have 
(u', V — u) + aoiu, V — u) + ai(u, u,v — u) + j{vr) — j{ur) — {f,v — u) 




— CTrVr)ds > 0, 



for all V G Vn,a- Hence m is a solution of Problem VIo-. 

Next, let M be a solution of Problem VIo-- Taking m ± u as a test function in 
p.ip . with arbitrary v € V^, wc find that 

{u ,v) + ao{u,v) + ai{u,u,v) ^ {f,v) (Vu G K). (3.3) 

By a standard theory (see Propositions 1. 1.1 and 1.1.2]), there exists unique 
peQ such that ([^ holds. Therefore, cr^ G {H^/^{TiYy is well-defined, and thus 

(u', v) + ao{u, v) + ai{u, u, v) + b{v,p) - {ar, Vr) ^1/2(^^1)'' = (/' (^"^ ^ ^«)- 
Combining this equation with p.ip . we obtain 

- ((Tr,Wr - 1tr)^l/2(ri)d < j 9{\Vt\ - \Ur\)ds (VuGl4,o-), (3.4) 

and as a result of triangle inequality, | (iTr , Wr)jLfi/2(p^)d | < /p^ g\vr \ ds for v G Vn,a- 
In view of Lemma ET^l ii). this implies that for rj G H^^'^iTi)'^ 

I {(^r,V)m/2(T^y I = I (^r,?7r)ffl/2(p^)d | < 1 1 ?7r 1 1 (Fi < 1 1 J?! I Li (Fi )<i • 

By a density argument, we can extend ar to an element of (L^(r)'')' such that 
I {'^r,v)Li(T,Y I ^ M\lI[t,y (Vt? G L\[T,Y). 

Since (^^(ri)'')' = ^^^^(ri)'*, we conclude |cr^| < g. Then ct^ • -f g\ur\ = 
follows from (j3.4l) with i; = 0. Hence {u,p) is a solution of Problem PDE. □ 
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3.2. Main theorem. Proof of uniqueness. We are now in a position to state 
our main theorem. We assume: 

(51) f &H\Q,T-L\nY). 

(52) g e ffi(0,T;L2(ri)) with 5(0) € H\Ti). 

(53) Mo e H'^i^Y n Vn.a, and SBCF is satisfied at t 0, namely, 

|o-r(uo)| < 5(0) and (Tr (uo)wor + 5'(0)|uor| = a.e. on Ti. 
Note that (Tr(uo) can be defined in a usual sense. 

Theorem 3.2. Under (9I1)~(3S1)) when d — 2 there exists a unique solution u of 
Problem W„ such that 

u e L°°(0, T; 14,.), u' e L°°(0, T; L''{nf) D L''{0, T; K,.). 

When d = 3, the same conclusion holds on some smaller time interval (0,T'). 

We call the solution in the above theorem a strong solution of Problem Via- 
First we prove the uniqueness of a strong solution. The existence will be proved in 
Section [3.41 after some additional preparations. 

Proposition 3.1. Ifui and U2 are strong solutions of Problem VIo-, then ui — U2. 

Proof. Taking v — U2 and v — ui 'm (j3.ip for ui and that for U2 respectively, and 
adding the resulting two inequalities, for a.e. t e (0, T) we obtain 

(u'l - u'2, ui ~ U2) + ao(ui - U2, ui - U2) 

< ai{ui,Ui,U2 - Ml) + ai{u2,U2,Ui - U2) 

= -ai(ui - U2,U2,ui - U2) - ai{u2,ui - U2,ui - U2). (3.5) 
We deduce from (|2.3p . together with Young's inequality, that 

|ai(ui - U2,U2,Ui - U2)\ < C\\ui - U2\\]^2t^^^y\\ui - '«2||^f(s2)d|l"2|li/i(f2)rf 

< - W2|lHi(n)<i + C'II"2|Ih1(0)^I|w1 - "2|li2(o)d, 
|ai(-U2,Ml - U2,Ui - U2)\ < C\\u2\\m{n)A\ui - U2||^f(f2)d||ui - W2||^(^(f2)d 

< - u2\\H^(^ny + C|I"2|Ihi(o)''II"i - Mh^ny- 
Combining (|2.ip and these estimates with p.Sp . we have 

^||"i - "2|li2(o)d < C'(||u2||^i(s;2)<i + \\u2\\Hi(ny)\\ui - U2\\l^ny- 
By Gronwall's inequality, we conclude 

ll"iW--2(0lli=(or <e^°''^"'"''-^-'^^''"^''-^<"'^^'*ll^ 

since mi(0) = 1*2(0) = uq. (Note that /q (Il'"2||^i(f2)d + ll'"2||^i(Q)d) dt remains finite 
because u e L°°{0,T; H\ny.) Thus ui{t) = U2{t). □ 

Remark 3.3. In the case of SBCF here, the last term of p. 51) vanishes, according to 
Lemma [2. 2 f i). We did not use that fact because we would like to make our proof 
of uniqueness remain unchanged when we deal with LBCF. 

Concerning the associated pressure, we find: 

Proposition 3.2. Under the assumptions of Theorem \'6.2[ let u be the strong so- 
lution of Problem VI^, and p be the associated pressure obtained in the proof of 
Theorem [3U Then p e L°°{Q,T;Q). 
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Proof. For a.e. t e (0,T), the well-known inf-sup condition (see [TUJ 1.(5.14)]), 
together with p.2p . (12. 4p . and |(Tt| < 5 a.e. on Fi, yields 

b{v,p) 



miL^n) < sup 



< 



lu'llL^^ny + C\\u\\Hi(ny + C||w||^i(jj)d + C\\g\\L2{ri) + \\f\\L^{n) 



Since RHS is bounded uniformly in t, p is in L°°(0, T;Q). □ 

3.3. Regularized problem. To prove the solvability of Problem VIo-, we consider 
a regularized problem VI^ , which is shown to be equivalent to a variational equation 
problem, denoted by Problem VE^. 

Before stating those problems in detail, for fixed e > we introduce 



where is a regularization of | • | having the following properties: 

(a) pe e C^(R'*) is a nonnegative convex function. 

(b) For all z e R'', it holds that 

\pe{z) - \z\\ < e. (3.6) 

(c) If tte denotes Vpe, for all z e M'' it holds that 

\a{z)\ < 1 and a^iz) ■ z>0. (3.7) 

(d) Let Pe denote the Hessian of p^, namely, /3e,ij — q^.q\. for j,j — 
Then /3e is semi-positive definite, that is, 

*y/3.(%>0 (V2/,zeM''), (3.8) 

where *y means the transpose of y. This is a consequence of the convexity of 

Such does exist; for example, let p^ be given by Pt{z) — \z\ — (1 — f )e if l^l > e, 
P(;(z) = ^(1 — cos ^|z|) if |z| < e. Then some elementary computation shows that 
Pe enjoys all of (a)-(d) above. 

Remark 3.4. One could use the Moreau-Yoshida approximation of | • | as pe, which 
is considered in [17], but it is only in C^K''), not in C''^{R'^). 

Since pe is differentiable, the functional je is Gateaux differentiable, with its 
derivative Dj^{rj) € (iJ^/^(Fi)'')' computed by 

je{r] + h£.) - jeiv) 



{DMv),Om^Hr,y - 1™ ' ^ " ^ / gaM ■ ^ds (3.9) 



for ?7,^ e i/i/2(Fi)''. 

We are ready to state the regularized problems mentioned above. 

Problem VI^-SBCF. For a.e. t £ (0,T), find u,{t) e Vn,a such that <(t) e 
L^{n)'^, UeiO) = 4 and 

(W^,U - Me) -I- ao(Ue,W - Me) + aiiUe,Ue,V- U^) + je(l'r) - jt{Ucr) 
>(/,«-Ue) (V«eK.a). (3.10) 

Problem VE^-SBCF. For a.e. t e (0,T), find u^it) g K,<t such that u'^{t) e 
L2(17)'^, m,(0) = and 

(u^,t;)-|-ao(u<:,w)-l-ai(we,We,w)-)- / gae{uer)-Vrds = (/, u) (Vu G K,o-)- (3.11) 
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Here, Uq is a perturbation of the original initial velocity uq. The way one obtains 
Uq from Mo is described later. By an elementary observation (e.g. [5J Section 3.3] 
or [27l Lemma 3.3]), we see that: 

Proposition 3.3. Problems VI^ and VE^ are equivalent. 

Now we focus on the construction of a perturbed initial velocity Uq. Since uq S 
H'^{QY satisfies SBCF by {3^, it follows from the Green formula ao{uo,v) = 
{—lyAuo^v) + /p^ (Jr{uo) ■ Vrds, for v £ Vn,a-, that 



ao{uo,v - uq) + J g{0)\vr\ds-J g{0)\uor\ds > {-vAuo.v - uq) 

{yveVn,.). (3.12) 

Here we consider the regularized problem: find Uq € Vn,a such that 



aoiuQ,V - Uq) + g{0)peiVr)ds ~ g{0)peiuQ^)ds > {-vAuq,V - Uq) 

(Vi^eK,.), (3.13) 

which is equivalent to (cf. Proposition [ 



ao{u'Q,v)+ 5(0)a,(u^^) -w^ds = (-i/Auo,^) (Vu G K..)- (3.14) 

"'El 

By a standard theory of elliptic variational inequalities , p.l3p admits a unique 
solution Uq, which is the perturbation of uq in question. With this setting, we find: 

Lemma 3.1. (i) When e — > 0, ^ strongly in H^{fl)'^. 
(ii) 4 e H^in)"^ and 

\\<\\HHny < CiWi^AuQW + ||<?(0)||j^i(r,)). (3.15) 

Proof, (i) Taking t; = uq in (|3.13p and v = Uq in (|3.12p . adding the resulting two 
inequalities, applying Korn's inequality, and using (13.61) . we conclude 



"Iko - woIIhuo)^ < / giO){\u'Q\- pe{u'Q))ds+ g{0){pe{uQ) -\uQ\)ds 
<2e g{0)ds->-0 (e^O). 

"'El 

(ii) Since g{0) € H^iTi) by (321)j we can directly apply the regularity result [?71 
Lemma 5.2] to the elliptic variational inequality p.l3p . and obtain (I3.15p . Though 
our pc and are different from those of |27j , it makes no difference in the proof of 
that lemma. □ 



Remark 3.5. (i) As a result of (i) above, for sufficiently small e > we have 

l|wol!L2(f2)<' < 2||uo|!L2(f2)<' and \\uQ\\Hi{n)i < ^uqWh^u)"- (3-16) 

(ii) Concerning the regularity of the domain, [27j assumes that Tq and Fi are 
class of and C* respectively, which is sufficient for our theory as well. 

Remark 3.6. In [27 , dealing with the stationary problem, the author stated that 
g g H^/^(ri) was enough to derive u € H'^{D,)'^ and p G H^{fl). However, it 
turned out that his proof presented there worked only for g e H^{Ti); see the 
errata by the same author. This is why we have assumed 5(0) e H^(Ti) in 
not g{0) e H^/^iTi). 
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3.4. Proof of existence. Due to Proposition l3.31 we concentrate on solving Prob- 
lem VE^. In doing so, we construct approximate solutions by Galerkin's method. 
Since Vn^a C H^{U)'^ is separable, there exist members wi,W2,-.- G Ki,ct, linear 
independent to each other, such that lJm=i span{«^fe}i!Li <^ ^n.a dense in H^{nY. 
Here e is fixed, and thus we may assume wi = Uq. 

Problem VE^"'-SBCF. Find G C^{[0,T]) (fc = 1, ...,m) such that Um & 
defined by Um — X^IILi Ck{t)wk satisfies u„i{0) — Uq and 



(fc = l,...,TO). (3.17) 

Since G C^iW^Y, the system of ordinal differential equations p. 171) admits 
unique solutions Ck G C^([0,T'])(/c = l,...,rn) for some T < T. The a priori 
estimate below shows T can be taken as T, so that we write T instead of T from 
the beginning. 

Proposition 3.4. Let (3T|)-(33]) be valid and e be small enough so that (|3.16p holds. 

(i) When d^ 2, u„, G L°-(0,r;K,.) and G L°^iO,T; L\n)'')nL^O,T;Vn,a) 
are bounded independently of m and e. 

(ii) When d — 3, the same conclusion holds for some smaller interval {0,T'), 
which can be taken independently of m and e. 

Proof. Due to space limitations, we simply write II m|| £,2, ||g||L2, H/Hl^, . . . instead 
of ||M||L2(o)d, \\g\\L^{ r,), U WL^n)", ■ ■ ■ and so on. 

(i) Multiplying (j3.17p by Cfc(i), and adding the resulting equations for k = 
1, m, we obtain 

(Kn^'^m) + aoiUm^Urn) + / gCteiUmr) ' "mrrfs = (/, Mm), 

where we have used Lemma I2.2f i) . It follows from (12. ip and (13. 7p that 

l^hrnWh +a||u„,||^i < (/,lim) < 1 1 / 1 1 ^2 || M„ 1 1 < ^||u,„||^l + ^^-11/11^2, 

2 at 2 2a 

which gives 

Consequently, for < t < T 



^^\\u^\\h+a\\u„,\\lr<C\\f\\l.. (3.18) 



2 , „, / IL, ||2 1, ^ |l„,e||2 , ri II /l|2 



+ « / WuAdt < holli2 + C / ||/||i2dt. (3.19) 

JQ JQ 

Since ||MollL2(f2)d < 2||uo||l2(s"2)<' by assumption, wc find that ||um||L°c(o,T:L2) and 
ll'"m||L2(o,T;V„ „) are bounded by C(/, uq) independently of m and e. 

Next, we differentiate (I3.17P with respect to t, which is possible because Ck{t)'s 
are in C^{[0,T]), to deduce 

{u['^,Wk) + ao(u^, Wfe) + ai(u^,u„, Wfe) + ai(u„, Wfe) 



+ / g' a^{umr) ■ Wkrds + I g^u'^^P^Wkr ds = {f ,vuk) (A: = 1, ...,m). 
Multiplying this by c'^.(t), and adding the resulting equations, we obtain 

("m: "In) + ao(u'„, <„) + Ol mr 

"'ri 

+ / g'u'^,l3,iu^r)u'„-,,ds^if,u'J, (3.20) 



ri 
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where we have again used Lemma I^T^ i). Here, 

) < C||w^„||i2||-u,„||^i||ii',J|^i (by (123)) 

< ^\WJ\m+C\\u,JUWJ\h, (3.21) 

<\\9'\\L4<r\\LHr,y (by (EH)) 

< C||.9'|U2||<J|^i (by Lemme.^i)) 
<^\WJ\l.+C\\g'\\l., 



/ g 'u'^^l3,{urnr)u'„,^ ds > 0, (by g > and ^) 



l(/',OI < Wf'h^u'JlH^ < + cWfWh. 



a 
6 

Collecting these estimates, it follows from p.20p that for < i < T 

|||<Jli.+a||<Jp^. <C7(||/'||i. + ||.g'||i.) + C||u™||^.||<Jli.. (3.22) 
If the second term of LHS is neglected, Gronwall's inequality leads to 

wKnimh < (^Hnmh +cj\\\n\h + \\9'\\h)dt^ e^^o n-^f.i'^*. (3.23) 

Provided that ||Wm(0)||^2 is bounded independently of m and e, estimate p.23p 
gives the boundedness of ||Mmlli°°(o, T;L2) because we know that of ||um||L2(o,r;V„ „) 
due to (|XT9)) . Then, by (pTS]) and ((XT9t we have 

(^\\um{t)\\]ji < C||/||^2 - |KJ|l2||w„,||l2 < C(/,g,uo), 

which implies ||w„j||/^oo(q y.y^ ^-j is bounded. Finally, integrating p.22p . we see that 
ll"JnlU2(o,T;y„,„) is also bounded. 

To show the boundedness of ||u^(0)|j^2, we multiply p.l7p by cj,.(t), add the 
resulting equations, and make t = 0, arriving at 



aoK>"m(0))+ / 5(0)a.«)-MLr(0)ds 

•/Ti 



= if{0),u'jO)). (3.24) 
From the construction of Uq, especially p.l4p . we have 

= |(-^.Auo,<„(0))| 

<C\\uo\\h4<M\\l-- (3.25) 
Furthermore, by Schwarz's inequality, Sobolev's inequality and p.l5L 

|ai«,uS,<„(0))| < C\\ulU^\\ul\\H4<nm\L^- < C\K\\U\<nmL^ 

<c{\\uo\\H^ + \\gmH^nuLm\L^- 

Combining these estimates with p.24p . we obtain 

||«™(0)|U2 < ||/(0)|U2 + C\\uo\\h^ + C{\\uo\\h2 + hmimf, 
which proves the boundedness of (0)11^2. This completes the proof of (i). 
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(ii) The discussion before p. 211) and the observation for ||wJ^j(0)||/,2 arc the same 
as (i). What changes from the case d = 2 is that when d = 3, instead of p.2ip . we 
only have (by (|2.3p and Young's inequahty) 

< C*! I ^^m II L^^ 1 1 1 1 ffi IKJ I ^/i^ 

< -/\\u„^\\h^\\u'J\hi+C\\u„,\\hi\\u',J12, 

for a constant 7 > which can be arbitrarily small. We choose 7 satisfying 
7l|Mo||fl-i < and from p.l6p we obtain 7||wq||j/i < j^. Let T' > 0, which may 
depend on to, e at this stage, be the maximum value oit such that 7||7im(t)||/ji < ^. 
If 7||u™(t)||^i < f for aU < t < T, we set T' = T. Since 7|ju„(0)||^i < f and 
Um(t) is continuous with respect to i, such T' does exist, and furthermore if T' < T 
then 7||7i„(i)||/fi = |. 

Therefore, in place of p.22p we obtain 

|||<Jli. + aWu'Jll, < C(||/'||i. + ||.g'||i.) + C\MH4<n\\h (0 < t < T'), 
which leads to the boundedness of ||Wm||L2(o.T';y„ and ||'Um||L°°(o,T';L2), together 

with ||'Um||L~(0,T';V„,„)- 

Finally, let us prove that T' is bounded from below independently of to and e. 
In fact, if T' < T then we see that 

-Z- < ||w„j(r')||^l - ||u„(0)||^fl < ||'U„i(T') - Ura{0)\\H^ = 
127 

< / W^UmH^dt < VT\\u'J\l2^o,T':V„^)- 

Jo 

Since we already know ||Wml|L2(-o y/.y^^ ^-j is bounded, we obtain the lower bound for 
T'. This completes the proof of Proposition 13.41 □ 

Remark 3.7. (i) A naive computation gives, by p.7p . 

< ll.9(0)||L2(ri)||w'„^(0)|li2(ri)d, 

but ||Mj„^(0)||i2(rj)£i cannot be bounded by ||u'„(0)||i2(Q')d in general. Therefore, 
the perturbation of uq, which is based on the compatibility condition in (331) j is 
essential in deriving (I3.25p . 

(ii) If 0? = 3 and /, g, uq are sufhciently small, we can prove j\\urn{t)\ \ h'^ (n)'' 

< ^ for all < i < T, and consequently the existence of a global solution. 

As a final step for our proof of the existence, we discuss passing to the limits 
TO — >■ (X) and e ^ 0. The proof below is valid for both d = 2, 3, except that when 
d = 3 we have to replace T with T' given in Proposition 13.41 

Proposition 3.5. (i) Under the assumptions of Proposition 13.41 there exists a 
solution of Problem VI^ such that all of \\u^\\i^oc(^Q T;Vn a)' ll''^ellL'^(o,T;V'„ „)j o-i^d 
ll'"ellL°°(o.T;L2(s7)<') bounded independently of e. 
(ii) There exists a strong solution of Problem VIo-. 

Proof, (i) As a consequence of Proposition 13.41 there exists some and a subse- 
quence of {um}$^=ii denoted again by {um}m=i, such that Ue G L°°{0,T;Vn,a), 

< e L^{0,T;Vn,a) f] L°°iO,T; L^ifl)"^), and when to -> 00 

u„i weakly-* in L°°(0, T; K,ct), 

u'^ < weakly in L^{0, T; Vn,a) and weakly-* in L°°(0, T; L'^{nf). 




u',nit)dt 



f giO)a,{ut,,)-u[^,{0)ds 
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We also find that all of ||ue||L~(o,T;y„,„), \\K\\L^{o,T:y„,^), and ||uellL~(o,T;L2(f2)i) 
are bounded independently of e. 

Let us prove solves Problem VI^- By Proposition 13. 31 it suffices to show that 
Ue is a solution of Problem VE^. Multiplying ([3T7)) by an arbitrary (j) £ C^{0,T) 
and integrating over (0,T), we obtain 

4>it)Uu'^,Wk) + aQ{u.,ri,Wk) + ai{Um,Um,Wk) + / gaeiujnr) ■ Wkrds 

dt^O {k ^ l,...,m). (3.26) 



(Lwk)} 



It follows from [29l Theorem III. 2.1] that the embedding 

{v\ve l2(0, T; H\n)'^), v' e L^{n X (0, T))"^} ^ l2(0, T; L^fl)'^) 

is compact, so that Um strongly in L^{Q^T\ L'^{^Y). Moreover, since the 

trace operator H^[n x (0,r)) — > L'^{Ti x (0,r)) is compact, u„„- — ?> u^r strongly 
in L^(ri x (0,T))'*. In particular, u^t Wer a.e. on Fi x (0,T), and thus the 
continuity of ae(z) yields a^{umT) — > cte(MeT) a.e. on Fi x (0,r). Making m — >■ cxd 
in (j3.26p . together with Lebesgue's convergence theorem, we see that 



0(i)|(u^,Wfe) + ao(u<:,Wfe) + ai(u£,U£,Wfe) + / ga^{u^r) ■ Wkrds 

-{f,Wk)}dt^O (fc = l,2,...). 

Since lJm=i span{u;fc}^j^ — Vn,a-, the above equation is valid for all test functions 
V S Vn^a- Hence p. lip holds for a.e. t g (0,T), which implies that is a solution 
of Problem VE^. 

(ii) As a result of (i), there exists some u and a sequence ei {I ^ oo), 
to which we drop the subscript I for simplicity, such that u e L°°{0,T;Vn^a), 
u' e L^{0,T;Vn^a) n L°°{0,T; L'^iiiy), and when e 

^ M weakly-* in L°°{0,T; 14, o-), 
< ^ u' weakly in L^{0,T;Vn,a) and weakly-* in L°° {0,T; L^{nf). 

As before, one sees that ^ u strongly in L^{0, T; i'*(i7)'^) and u^r — >■ strongly 
in L^(Fi X (0, T)). In addition, ^ u weakly in i^(0, T; V^^o-), and thus it follows 
that /J" ao(w, u) < lim ^ ,q JJ" ao(we, We) di. 

Following the technique of p. 56], we let v{t) £ L^(0, T; 14,^) be arbitrary. For 
a.e. t e (0,T), we take v = v{t) in p.lOp and integrate the resulting equation over 
(0, T) to deduce 

W — -Ue) -I- ao(Me, V — U^) + ai{u^, Me, V — Ue) 

+ je{ir) - jeiUer) "(/,«- We)}dt > 0. (3.27) 

In view of p.6p . together with triangle inequality and Lipschitz continuity of j, 

rj-\ rj^ rj-t rj-i 

we have je{vT)dt — j{vr)dt and /p je{ugr)dt — ^ j{ur)dt when e — > 0. 
Therefore, taking the lower limit lim ^. ,q in p.27p gives 

^{u' ,v — u) + ao(u, V — u) + ai(u, m, u — u) -f- jX^t) ~ jiur) ~ {f,v — u)^dt > 0. 

A technique using the Lebesgue differentiation theorem allows us to conclude that 
u satisfies p.ip at a.e. t — to (for more detail, see 5, p.57]). 
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Concerning the initial condition, since the trace operator H^{fl x (0,T)) — > 
L^{nx{0}) is continuous, Lemma lXlT i') leads to u(0) = linie_s.o Ue(0) = linie^.o = 
mq. Hence u is a strong solution of Problem VIo-. □ 

Propositions 13.11 and I3.5r ii') complete the proof of Theorem 13.21 

4. Navier-Stokes Problem with LBCF 

4.1. Weak formulations. Throughout this section, we assume / e L^{ilx (0, T)), 
uq G Vr.a, and g G L'^iTi x {0,T)). Further regularity assumptions on these data 
will be given before Theorem 14. 2 1 In addition, the barrier term j„ is simply written 
as j. A primal weak formulation of (ll.ip - (ll.4p with (jl.6p is as follows: 

Problem PDE-LBCF. For a.e. t e (0,r), find {u{t),p{t)) e Vr x Q such that 
u'{t) e L'^{Q)'^, u{0) = uq, (t„ is well-defined in the sense of Definition l2.11 |cr„| < g 
a.e. on Fi, and (TnUn + g\un\ = a.e. on Fi. 

Remark 4.1. More precisely, "|cr„| < g" implies that an G {H^^'^{Ti)y actually 
belongs to ^^^(ri) with ||cr„||ioo (r^) < 1. In particular, cr„ G L^(ri). 

Throughout this section, we refer to Problem PDE-LBCF just as Problem PDF. 
Similar abbreviation will be made for other problems. Next, as in SBCF, we pro- 
pose: 

Problem VI-LBCF. For a.e. t e (0,r), find {u{t),p{t)) e Vr x Q, such that 
u'{t) e L^{n)'^, u{0) = ua and 

(w', V — u) + ao{u, V — u) + ai(u, u,v — u) + b{v — u,p) 

+j{Vn) - j{Un) >if,V~ U) (Vv G K), (4.1) 

. biu,q)^0 (Vqeg). (4.2) 

Problem VI^-LBCF. For a.e. t G (0,r), find u{t) e Vr,a such that u'{t) e 
L^{n)'^, u{0) = uo and 

(m', W - m) -I- ao{u, V-u)+ ai{u, U,V -u) +j{Vn) - > (/; - ") (Vw e T^r,o-)- 

(4.3) 

Unlike the case of SBCF, Problem VIo- is not exactly equivalent to Problem 
PDE, as is shown in the following theorem. 

Theorem 4.1. (i) If {u,p) solves Problem PDE, then u solves Problem Yl^. 

(ii) Ifu solves Problem VIo-, then there exists at least one p such that {u,p) solves 
Problem PDE. // another p* satisfies the same condition, then for a.e. t S (0,T) 
there exists a unique <5(i) S M such that 

p{t)^p*{t)+6{t) and an{u{t),p{t)) = an{u{t)y-{t))-5{t). (4.4) 

(iii) In (ii), if we assume furthermore Un(t) ^ 0, then S{t) = 0. Namely, the 
associated pressure is uniquely determined. 

Proof, (i) This can be proved by the same way as Theorem 13. II 

(ii) For a.e. t e (0,r) and v e it follows from that {u',v) + ao{u,v) + 
ai(u, u, v) — (/, v), and thus there exists unique p € Q such that 

{u\v) + ao(u,u) + ai{u,u,v) + b{v,p) = {f,v) (Vw E V). 

According to Definition 12.11 (7„ = cr„(u,p) is well-defined, so that 

{u',v) + ao{u,v) + b{v,p) + ai{u,u,v) - {cr„,Vn) ^1/21^^1) = (/' ^) C^^ ^ K)- 
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Substituting this equation into (|4.3p . we obtain — ((7„,u„ — Un) H'-/^{ri) — ji'^n) — 
j{un) for all V E Vr^a- It follows from Lemma [2.4^ 11') that 

\{^n,v)m/2^r^)\< f 9\v\ds {Vrj e h'^'{t,) n Ll{r,)). 

The Hahn-Banach theorem allows us to extend (T„ to a linear functional cr„ : 
Lg(ri) — > R satisfying the same inequality as above for all G Lg{Ti). There- 
fore, (T„ e L'^^g{Ti) and |cr„| < g. In addition, (T„u„ + (7|u„| = follows. 

Since (t„ — (t„ vanishes on H^^^{Ti)nLQ{Ti), there exists a constant S{t) such that 
o'n — <7n — <5(i). Now, by setting p{t) — p{t) + 6{t), it follows that (T„ given above 
actually equals an{u{t) , p{t)) and that {u{t),p{t)) solves Problem PDE. Relation 
(|4.4p can be verified by a similar argument. 

(iii) Since /p^ u„ ds = J^divudx = 0, the assumption u„(i) ^ implies that 
there exist subsets A+, A_ of Fi with positive d—1 dimensional Lebesgue measure 
satisfying Un{t) > on and < on A^. Because |cr„| < g and (Jn^tn + 

g\un\ = on Fi, (t„ — —g{t) on ^+ and cr„ = ^(t) on Hence 5{t) in (j4.4p 
cannot be other than zero. □ 

Remark 4.2. Since |cr„| < (5(i) is no more than 2g{t) nor less than —2g{t). 

4.2. Main theorem. Let us state our main theorems for the case of LBCF. As in 
SBCF, some compatibility condition is necessary; it is rather complicated because 
normal stress at t = involves a pressure at i = 0, which is not given as a data. 
The precise description is as follows: we say that LBCF is satisfied at t = if 
Mo G H^{fl) n Vr,a and there exists po e H^{fl)'^ such that 

|o-„(-uo,_po)| < g(0) and cr„(Mo,_po)"Ori + 5(0)|uori| = a.e. on Fi (4.5) 

We remark that a similar compatibility condition appears in nonlinear semigroup 
approaches (see [ZlIB])- 

Furthermore, in order to overcome a difficulty arising from (|1.7p . we need no- leak 
condition at i = 0, that is, uqu = on Fi. In view of (j4.5p . this is automatically 
satisfied if |cr„(uojPo)| < 5(0) on Fi. Examining our proof of the a priori esti- 
mates carefully, one finds that this assumption can be weaken to the condition that 
ll?^0n||L2(ri) is sufficiently small. 

Including what we have discussed above, we assume the foUowings: 

(LI) f eH\0,T;LHn)''). 

(L2) g G H^{0,T;L^{Ti)) with 5(0) G H^Ti). 

(L3) uo G H'^in)'^ n Vr^a, and LBCF is satisfied at t 0. 

(L4) UQn — a.e. on Fi. 

Theorem 4.2. Under (L[T|)-(I|31) above, there exists a unique solution u oj Problem 
VIcr on some interval (0,T'), with T' < T, such that 

u G L°°(0, T'; K,a), u' G L°°(0, T' ■ L^{9.f) n £^(0, T'; Vr^o). 

The uniqueness can be proved by the same way as Proposition 13.11 We can 
also obtain p G L°°(0, T'; L^(ri)) by a similar manner to Proposition [321 using the 
rather infamous inf-sup condition (see [271 Lemma 2.2]) 

CbllL^(o) < sup TT^P^ (Vp G L\n)). 

The rest of this section is devoted to the proof of the existence. To state regu- 
larized problems, for fixed e > we introduce 

jciv) = / gpe{v)ds, 
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where is a function such that 

(a) pe € C^(R) is a nonnegative convex function. 

(b) For ah z e R, |p,(z) - |z|| < e. 

(c) If Qfe denotes dp^jdz, for all 2; G M it holds that 

|q;(z)| < 1 and ae(z)z > 0. (4.6) 

(d) Let /3e = d'^p^/dz^. Then /3e > (due to the convexity of p^). 

Such does exist; for example, if we define Pt{z) — \z\ — (1 — |-)e if \z\ > e, 
Pt{z) = ^(1 — cos ■^\z\) if |z| < e, then this pe enjoys all of (a)-(d) above. 

Since Pf is differentiable, the functional je is Gateaux differentiable, with its 
derivative Dj^{r]) G H^^^lTi)' computed by 

{D3M,Om/^ro^ f a^M^ds (V77,eGi?'/'(ri)). 

Now let us state the regularized problems. 
Problem VI^-LBCF. For a.e. t e (0,T), find Ue{t) e K.<t such that u'^{t) G 
L^{n)'^, UeiO) = 4 and 

(m^,U - + ao(u£,W - Ue) + ai{u^,Ue,V- Me) + ^^(Wn) - jeiUen) 

>(/,«-«.) (Vi;eK,a). 

Problem VE^-LBCF. For a.e. t € (0,T), find u^{t) € Vr,a such that <(t) e 
L2(17)'^, m,(0) = 4 and 

{u'^,v) + ao{u^,v) + ai{u^,u^,v) + ga<:(Me„)u„ds = (/, w) (Vu G K,<t). 



As in Proposition 13.31 Problems VI^ and VE^ are equivalent. The construction 
of the perturbed initial velocity Uq is similar to that of SBCF. In fact, since LBCF 
holds at t = by (1(3]), the Green formula leads to 

ao{uo,v-uo)+ g{0)\vn\ds-- g{0)\uon\ds 

We consider the regularized problem; find Uq G V-r,a such that 

ao(uo,w-Uo)+ J 9{0)Pe{vn)ds- J g{Q)pe{uQ^)ds 
>{~uAuo+Vpa,v-ul) {yveVr,a), (4.7) 



which is equivalent to (cf. Proposition 13.31 



ao(wo, v) + / g(0)a,«„)w„ ds = {-i^Auq + Vpo, v) (Vw G Vr,a)- (4.8) 

The elliptic variational inequality (|4.7p admits a unique solution Uq, which is the 
perturbation of uq in question. With this setting, we find; 

Lemma 4.1. (i) When e 0, — > mo strongly in H^{D,)'^. In particular, it 
follows that Uq in i^(ri). 
(ii) ul G H^{n)'^ and 

hollff^(o)- < C{\\vAuo + "^PoWmny + \\giO)\\m{r,))- (4.9) 

Proof, (i) is proved by the same way as Lemma [3.1f i). Since g{0) G H^{ri) by 
(I®, (ii) is a direct consequence of [271 Lemma 4.1]. □ 
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Remark 4.3. By (i) and (L|l]), for sufficiently small e > we have 

(4.10) 

where a and 71 are the constants in (|2.ip and (|2.5p respectively. 



Remark 4.4. As in SBCF, if Tq is and Fi is C^, then we can apply Lemma 4.1 
of [37]. On the other hand, 5(0) G H^^^(Ti), stated in ^7\, is actually insufficient 
to deduce the H'^-H^ regularity (see the errata of [27]). 

To solve Problem VE^, let us construct approximate solutions by Galerkin's 
method. Since VV,(j C H^{il,)'^ is separable, there exist wi,W2,... € Vr.a: linear 
independent to each other, such that lJm=i spanjwfelJJ^^ C Vr^a dense in H^{n)'^. 
Here we may assume wi = Uq. 

Problem VE^"^-LBCF. Find Ck e C^i[0,T]) {k ^ 1, ...,m) such that Um e Vr,a 
defined by Um — X^fcLi Ck{t)wk satisfies Wm(0) = Up ^^"^ 

(-U^,Wfc) +ao(Um,Wfc) +ai(M„i,M™,Wfc) + j gae{Umn)wknds={f,Wk) 

{k = l,...,m). (4.11) 

Since G C^(M), there exist unique solutions Ck € C'^([0,T])(fc = l,...,m) for 
some T, which may depend on m and e at this stage. 

Proposition 4.1. Assume (L[I])^(IJll)i cind let e > be sufficiently small so that 
(j4.10p holds. Then there exists some interval {0,T') such that Um G L°°(0,T';Vr^a) 
and mJ„ G L°°(0, T'; L^(r2)'^) H L^(0, T'; V^^o-) are uniformly bounded with respect to 
m and e. Here, T' is independent of m and e. 

Proof. Due to space limitations, we sometimes simply write ||u||x,2, ||(7||i2, in- 
stead of ||w||i2(s^)2, ||g||L2(ri), when there is no fear of confusion. 

First we consider the case d — 2. Multiplying (14. lip by Ck{t) for k — 1^ ...,m, 
adding them, using (|2.1I) . (|2.5p and (|4.6|) . we obtain 

^^Il'^m||i2 + (a -7l||u„m||L2(ri))||u™||^l < if,u,n)- (4.12) 

Since ||um„(i)||^2(pj-) is continuous with respect to t and (j4.10p holds, there exists a 
maximum value Ti e (0,r] of t such that 7i||umn(^)||L2(ri) ^ If this inequality 
holds for aU < t < f, we take Ti = f. Noting |(/,u„)| < f ||u,„||^i + ^||/||i2, 
we find from (I4.12p that 

j^\\u,n\\l2 + a\\u„,\\l^ <C\\f\\l. (0<t<Ti). 

Hence Mm £i°°(0, Ti;L'^) n L^(0, Ti; K.a) is bounded independently of to, e. 

Next, differentiating (|4.1ip . multiplying the resulting equation by c'j,(i), and 
adding them, we obtain 

(Wm,Mm) + ao «„,<„) + ai(u^,u„,M^) + fli (u„ , , uj„) 

+ / g'a,{umn)u'„,,^ds+ gl3^{umn)\u'^n\'^ ds = {f',u'„J. (4.13) 
Here, we estimate each term in (|4.13p as follows: 

)| < C\\u'„^\\L2\\Um\\H^\\Kn\\L^ 

<^\\u'J\''m+C\\u^\\U\<n\\L^: (4.14) 



\ai{um,u'^,u'^)\ 



r. 



< 7l||"mn||L2(ri)||W™|lHi < -^\\u',n\\m, 
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'ri 



<C\\g'h4uL\\m<Y^\\<\\m+Cy\\h, 



g/3<:(u„„)|u^„P ds > 0, 

\{f',u'J\<^\\u'J\l,+C\\f\\l.. 

Collecting these estimates, we derive from (I4.13P that for < t < Ti 

|||<JU.+a||^:„|p^. <C'(||/1|i. + ||.9l|iO + C||u„||?,.||<J|i.. (4.15) 

Combining the technique used in Proposition 13.41 with (|4.8p and (|4.9p . we ob- 
serve that |KJ|i=e(o,Ti:L2), IK„IIl2(o,Ti;K.<,), and ||wm||L-(o,Ti;y,.^) are bounded 
by C{f,g,uo,Po)- 

It remains to show that Ti is bounded from below independently of to, e. If 
7i|l'"mTi(Ti)||L2(p^') < a/4 and thus Ti = T, we can extend Um{t) beyond t — T and 
repeat the above discussion until we reach either 

^max^7i||u™„(i)||i2(rj) < q:/4 or 7i||u™„(ri)||i2(rj) = a/4. 

In the former case Ti = T. In the latter case, we have 
a 

Z < Il"mn(ri)||i2(ri) - ||Umn(0)||L2(ri) < ||Umn(T'l) - "mn (0) || L2 (p ^ 

»7l 

< / \\u'mn{t)\\L^ri)dt < C \\u'^\\Hi{nydt < C^/J\\\u'J\L2(^o^T^■,V^,^)■ 

Hence Ti is bounded from below, and we complete the proof for d = 2. 

Second let us consider the case c? = 3. What changes from d = 2 is that (|4.14p 
is replaced with 

\ai{u'^,u^,u'^)\ < C\\u'^\\]^^^\\um\\HA\u'm\\m 

< l2\\Um\\m\\u'^\\Hi+C\\Um\\H^\\Kj\l2, 

where 72 can be arbitrarily small. We choose 72 satisfying 72||uo||ifi < so that 
72||wo||ffi < ^ by virtue (|4.10p . Let T2 be the maximum value of i e (0,T] such 
that 72||Mm(i)||//i < j^- If this inequality holds for aU t e (0,T], we set T2 = f. 
Such T2 does exist, and if T2 < T then -f2\\um{T2)\\ = j^. 
Therefore, setting T' = min(Ti,T2), instead of (|iTT5|) we get 

j^\\u'J\l^ + aWu'Jll, < C(||/'||i2 + ||5'|li2) + C\\uJ\H4<n\\h (0 < ^ < T'). 

As a consequence, we see that |i<„||L2(o,T';y,,„), \\u'„i\\l°°(o,T';L2), ||wm||L°-(o,T';y^.^) 
are bounded hy C{f, g,uo,Po)- 

Now, if Ti < T or T2 < T then T' are bounded from below as follows: 

a f^' 

< ||Wm«(r')||L2(ri) - l|Wm«(0)||L2(ri) < / I K„J I (Fi ) 



1271 

<C [ WKjH^dt < CVT\||<J|L^(o,T';y.,„) 
Jo 

i-T' 



2472 



< \\u^{T')\\m - WuUmm < I \\u'J\mdt < Vt\\u'J\l^o.t';V^.„). 

Jo 



When Ti — T and T2 = T, we can extend Um{t) beyond t = T and repeat the 
above discussion. This completes the proof of Proposition 23] □ 
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The last step of the proof, namely, passing to the limits m — >■ cx) and e — >■ can 
be carried out by the same way as Proposition 13.51 with n replaced by r and vice 
versa. This proves that a solution of Problem VIo- exists, which, combined with the 
uniqueness result, completes the proof of Theorem 14.21 

Remark 4.5. At first glance one may think Theorem 14. 2 1 where we get only a time- 
local solution in spite of a smallness assumption on uq even if d = 2, is too poor. 
However, in view of the fact that we obtain only time-local solutions in 2D case 
under the linear leak b.c. (see Theorem 6] or [H]), such limitations cannot be 
avoided to some extent. 

Remark 4.6. Under additional smallness assumptions on the data /, uo,poj we 
can derive global existence results for both d = 2 and d = 3. 

5. Concluding Remarks 

By the discussion presented above, we have established the existence and unique- 
ness, while we did not get in touch with higher regularity such as 

u e L°°{0,T;H\nf), p e L°°{0,T;H\n)). 

This is because some regularity results for the elliptic cases are not available. For 
instance, Problem VI(j-SBCF is rewritten as 

ao{u,v - u) + j{vr) - j{ur) > {f,v - u) - {u' , V - u) - ai{u,u,v- u) 

with F{t) e LP{il)'^ for some p < 2. If we prove this elliptic variational inequality 
has a unique solution in W'^'P{il)'^ when p < 2, then a technique similar to [251 
Theorems III. 3. 6 and III. 3. 8] allows us to deduce u{t) G H'^{ft)'^. Thereby, we need 
to extend the regularity theory of [22] to cases p^2. 

Acknowledgements 

The author would like to thank for Professor Norikazu Saito for encouraging him 
through valuable discussions. This work was supported by Grant-in- Aid for JSPS 
Fellows and CREST, JST. 

References 

[1] R. An, Y. Li, and K. Li, Solvability of Navier-Stokes equations with leak boundary conditions, 
Acta Math. Appl. Sin. EngL Ser. 25 (2009), 225-234. 

[2] M. Ayadi, M. K. Gdoura, and T. Sassi, Mixed formulation for Stokes problem with Tresca 
friction, C. R. Math. Acad. Sci. Paris 348 (2010), 1069-1072. 

[3] A. Yu. Chebotarev, Modeling of steady flows in a channel by Navier-Stokes variational in- 
equalities, J. Appl. Mech. Tech. Phys 44 (2003), 852-857. 

[4] L. Consiglieri, Existence for a class of non-newtonian fluids with a nonlocal friction boundary 
condition. Acta Math. Sin. (Engl. Ser.) 22 (2006), 523-534. 

[5] G. Duvaut and J. L. Lions, Les inequations en mechanique et en physique, Dunod, 1972. 

[6] H. Fujita, A mathematical analysis of motions of viscous incompressible fluid under leak or 
slip boundary conditions, RIMS Kokyilroku 888 (1994), 199-216. 

[7] , N on- stationary Stokes flows under leak boundary conditions of friction type, J. Corn- 
put. Math. 19 (2001), 1-8. 

[8] , A coherent analysis of Stokes flows under boundary conditions of friction type, 

J. Comput. AppL Math. 149 (2002), 57-69. 

[9] H. Fujita, H. Kawarada, and A. Sasamoto, Analytical and numerical approaches to stationary 
flow problems with leak and slip boundary conditions, Lecture Notes in Num. AppL AnaL 14 
(1995), 17-31. 

[10] V. Girault and P. A. Raviart, Finite element methods for Navier-Stokes equaitons, Springer- 
Verlag, 1986. 

[11] R. Glowinski, Numerical methods for nonlinear variational problems. Springer- Verlag, 1984. 



20 



TAKAHITO KASHIWABARA 



[12] J. G. Hcywood, R. Rannacher, and R. Turek, Artificial boundaries and flux and pressure con- 
ditions for the incompressible Navier- Stokes equations, Internat. J. Numer. Methods Fluids 
22 (1996), 325-352. 

[13] T. Kashiwabara, On a finite element approximation of the Stokes equations under boundary 
conditions of friction type. Part I: slip boundary problem, submitted (2011), 1-30. 

[14] , On a finite element approximation of the Stokes equations under boundary conditions 

of friction type. Part II: leak boundary problem, submitted (2011), 1-27. 

[15] H. Kawarada, H. Pujita, and H. Suito, Wave motion breaking upon the shore, GAKUTO 
Internat. Scr. Math. Sci. Appl. 11 (1998), 115 159. 

[16] N. Kikuchi and J. T. Odcn, Contact problems in elasticity, SIAM, Philadelphia, 1988. 

[17] D. S. Konovalova, Subdifferential boundary value problems for Navier-Stokes evolution equa- 
tions. Differ. Equ. 36 (2000), 878-885. 

[18] O. A. Ladyzhenskaya, The mathematical theory of viscous incompressible flow, Gordon and 
Breach, 1969. 

[19] C. Le Roux, Steady Stokes flows with threshold slip boundary conditions. Math. Models 

Methods Appl. Sci. 15 (2005), 1141-1168. 
[20] C. Le Roux and A. Tani, Steady solutions of the Navier-Stokes equations with threshold slip 

boundary conditions, Math. Mcth. Appl. Sci. 30 (2007), 595-624. 
[21] Y. Li and R. An, Semi-discrete stabilized finite element methods for Navier-Stokes equations 
with nonlinear slip boundary conditions based on regularization procedure, Numer. Math. 117 
(2011), 1-36. 

[22] , Two-level pressure projection finite element methods for Navier-Stokes equations 

with nonlinear slip boundary conditions, Appl. Numer. Math. 61 (2011), 285-297. 

[23] Y. Li and K. Li, Penalty finite element method for Stokes problem with nonlinear slip bound- 
ary conditions, Appl. Math. Comput. 204 (2008), 216-226. 

[24] , Locally stablized finite element method for Stokes problem with nonlinear slip bound- 
ary conditions, J. Comput. Math. 28 (2010), 826-836. 

[25] , Pressure projection stabilized finite element method for Navier-Stokes equations with 

nonlinear slip boundary conditions, Computing 87 (2010), 113-133. 

[26] S. Marusic, On the Navier-Stokes system witli pressure boundary condition, Ann. Univ. Fer- 
rara Scz. VII Sci. Mat. 53 (2007), 319-331. 

[27] N. Saito, On the Stokes equation with the leak or slip boundary conditions of friction type: 
regularity of solutions, Publ. Res. Inst. Math. Sci. 40 (2004), 345-383. 

[28] H. Suito and H. Kawarada, Numerical simulation of spilled oil by fictitious domain method, 
Japan J. Indust. Appl. Math. 21 (2004), 219-236. 

[29] R. Temam, Navier-Stokes equations, North-Holland, 1977. 

Graduate School of Mathematical Sciences, The university of Tokyo, 3-8-1 Komaba, 
Meguro, Tokyo 153-8914, Japan 

E-mail address: tkashiwaSms.u-tokyo.ac.jp 



